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ABSTRACT 

Context. Molecular hydrogen (H?) gas is the most abundant molecule in the circumstellar (CS) environments of young stars. It is thus 
a key element in our understanding of the evolution of pre-main sequence stars and their environments towards the main sequence. 
Aims. At the present time, little is known about the gas as compared to the dust in the environments of young stars. We thus observed 
molecular hydrogen around a sample of pre-main sequence stars in order to better characterize their gaseous CS environments. 
Methods. The FUSE (Far Ultraviolet Spectroscopic Explorer) spectral range offers access to electronic transitions of H2. We analyzed 
the FUSE spectra of a sample of Herbig Ae/Be stars (HAeBes) covering a broad spectral range (from F4 to B2), including the main- 
sequence A5 star/3 Pictoris. To better diagnose the origin of the detected molecular gas and its excitation conditions, we used a model 
of a photodissociation region. 

Results. Our analysis demonstrates that the excitation of Ht is clearly diiferent around most of the HAeBes compared to the interstellar 
medium. Moreover, the characteristics of H2 around Herbig Ae and Be stars give evidence for different excitation mechanisms. For 
the most massive stars of our sample (B8 to B2 type), the excitation diagrams are reproduced well by a model of photodissociation 
regions (PDR). Our results favor an interpretation in terms of large CS envelopes, remnants of the molecular clouds in which the 
stars were formed. On the other hand, the group of Ae stars (later than B9 type) known to possess disks is more inhomogeneous. In 
most cases, when CS H2 is detected, the lines of sight do not pass through the disks. The excitation conditions of Ht around Ae stars 
cannot be reproduced by PDR models and correspond to warm and/or hot excited media very close to the stars. In addition, no clear 
correlation has been found between the ages of the stars and the amount of circumstellar H2. Our results suggest structural differences 
between Herbig Ae and Be star environments. Herbig Be stars do evolve faster than Ae stars, and consequently, most Herbig Be stars 
are younger than Ae ones at the time we observe them. It is thus more likely to find remnants of their parent cloud around them. 

Key words, stars: circumstellar matter - stars: formation - stars: pre-main sequence - ISM: molecules 



1. Introduction 



Gradv et aLll200U lAugereau & Papaloizoull2004t iLagage eTaP 
2006h . Most evidence of disks around Herbig Ae/Be stars is 



A long standing problem m modern astrophysics is how stars j^^j^g^j. blueshifted absorption Unes, polarization, mass esti 
and planets form from their parent moleculai- clouds. It is gener- ^^^^^ f^^^ millimeter observations, extinction measurements, 
ally accepted that the collapse of an interstellar molecular cloud ^^^ ^- ^^ disk-model fitting to infrared specti-al energy distri- 
core to form a protostar naturally produces a disk-shaped neb- butio ns (IR SE Ds) (Bastien & Menard! 1 1990t iHiUenbrand et all 
ula in which planets might form. In the past twenty years, both g^g^. ^^^^^fan & Rav.. 19971) . Numerous observations have re- 
planets and disks have been observed around nearby stars. Some ^^^1^^ elongated structures around several HAeBes, with ve- 
youn^ stars with disks, e.g., the T Tauri star GM Aur ( |Riceetal. | j^^jjy gradients along the major axes characteristic of gas in 
120031), are suspected of harboring young planets, and recently a Keplerian rotation (Mannings & Sargent 1997, 2000). In addi- 
massive planet orbiting TW Hya, a young star surrounded by a ^-^^^ ^^e SEDs of some HAeBes show dips near lO^/m that can- 
well-studied disk, has beeii discovered by means radial velocity ^^^ ^e explained by normal dust-removal processes, such as the 
techniques (j Setiawan et alJ | 2008D . Poynting-Robertson drag and radiation pressure from the star It 
Disks have been directly iinagecl around T Tauri stars instead suggests that the dust structure is evolving because of the 
(e.g., HK Tau, [ Stapelfeldt et al. | |I998D, while direct or coro- ^^^y. breaki ng up. One p ossibility is that this is caused by planet 
nagraphic imaging of Herbig Ae/Be star disks has been per- formation (iBouwmanet al. 2003). In this context, HAeBes ai-e 
formed successfully only in a few cases (e.g.| Pantinet al.||2000; possible precursors of/? Pictoris and Vega-type stars, whose CS 

debris disks are believed to host planetary bodies. This raises 
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the possibility that the environment around the HAeBes truly 
represents a very early phase of planet formation. The physical 
parameters of HAeBes disks can thus be used to constrain the 
duration of the protoplanetary phase and the age at which the 
signatures of planets become visible in CS disks. 

On the other hand, mid-IR observations at high spatial reso- 
lution have revealed that the emission observed from HAeBes 
is generally not confined to optically thick disks but instead 
comes fr om more complex en vironments such as remnant en- 
velopes (Pol omski et alj|2002|). On the b asis of millimeter in- 
terferometric measurements, iNatta et all (2000) conclude that 
CS envelopes are much more common around Herbig Be stars 
(HBes) than around Herbig Ae/B9 stars (HAes). These con- 
clu sions have been con firmed by near-IR speckel interferome- 
try JLeinert et al.ll200l1) . This indicates structural differences be- 
tween Herbig Ae stars and Herbig Be stars. These results are 
fully compatible with the faster evolution of the more massive 
HBes. Indeed, due to their stronger radiation fields, the CS envi- 
ronments around HBes evolve faster than around HAes, which 
translates especially into a rapid depletion of the CS disk mate- 
rial. Consequently, the average HBe stars are younger than HAes 
at the time we observe them; therefore, it is more likely to find 
larger amounts of remnant CS material around Herbig Be stars, 
since less time is available to dissipate it. In this scenario, it is 
expected that the CS material is mostly concentrated at the outer 
edges of the close CS environment and distributed in more or 
less spherically envelopes. 

Molecular hydrogen, from which giant planets are primar- 
ily believed to form, is the most abundant molecule in the CS 
environments of young stars. The detection of H2 provides the 
most direct information about the gaseous content of the CS en- 
vironments of HAeBes, setting limits on the timescales for the 
dissipation of CS matter and possible planet building. Tentative 
detections of pure rotational H2 lines at infrared wavelengths 
towards a few Herbig stars and T Tauri stars h ave been pre- 
viously obtained with ISO-SWS dThi et al.ll20"ol . But ground- 
based observations have shown that contamination from sur- 
rounding cloud material can be important and that ISO coul d 
not have detected disk gas (" Richter et al.l2002l:[Sako et al.120051) . 
The H2 mid-IR emission lines have been clearly de tected in the 
disks of only two Her big stars, namely AB A ur JBitner et alj 
I2007h and HD 97048 dMartin-Zaidi et al.ll2007h . These studies 
show that the gas has not dissipated in the inner parts of these 
disks with ages of about 3 Myrs and that peculiar H2 excita- 
tion conditions are present in the disks. The electronic transi- 
tions of H2, which are much stronger than the rotational ones, 
occur in the far-ultraviolet (FUV) spectral range, a domain cov- 
ered by FUSE (Far Ultraviolet Spectrocopic Explorer). Previous 
FUSE obs ervations of selected you ng sta r s harboring CS disks 
dLecaveher des Etangs et al. 200ll 120031: iRobergeet ani200lh 
have demonstrated that such a study can make significant contri- 
butions clarifying issues raised in previous work. 

We present the results of our systematic analysis of the FUSE 
spectra of a large group of Herbig HAeBe stars with a variety of 
spectral types and Vega-type stars, including/? Pictoris. This pa- 
per is organized as follows. In Sect. |2] we present our sample 
stars and recall the methods used for data reduction and analy- 
sis, which have been presented in previous papers. The results 
obtained from the H2 analysis are given in Sect. [3] The detailed 
analysis and interpretation of our results for each group of stars 
are presented in Sects. |4]and|5] We explore possible evolutionary 
trends in our results in Sect. |6] Our conclusions are discussed in 
Sect. 12] 



2. The sample stars 

To perform this study, we have compiled a set of observations 
from different programs we carried out over the past six years, 
supplemented by some data from the archives (Table [T]). This 
results in a sample of 18 HAeBes with spectral types from F4 
to B2 and one main-sequence star, /? Pictoris, prototype of the 
Vega-type stars. Two stars, HD 141569A and HD 109573, be- 
long to the so-called transitional class and are passing from the 
pre-main sequence Herbig star stage to the zero-age mai n se- 
quence (ZAMS) stars. Both are known to possess CS disks dJural 
II99II; lAugereau et all 1 19991) with ages close to about 10 Myrs 
(TableH. 
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Fig. 1. The stars of our sample are plotted on this HR diagram. 
The evolutionary tracks between the pre-main sequence phase 
and the ZAMS are shown for th e Ae/B9 stars. For this plot, we 
used the Palla & Stable^ d 19931) pre -main sequ e nce tracks and 
the interpolation routines written bv lTesti et alj dl998.) . For the 
more massive Be stars, the pre-main sequence stage is less obvi- 
ous due to the faster evolution of these stars, which makes them 
appear directly on the ZAMS. 



Our sample also includes HD 135344. The spectral type of 
this star has not been determined precisely in the literature, es- 
timates ranging fr om F4 to AO, but it has been clas s ified as a 
Herbig star (e.g. The et al J [19941 iMalfait et all 1 19981: 1™!?^] 
2001). In the rest of this paper, we thus include this star in the 
sub sample called Herbig Ae/B9 stars. 

For several stars, some early results of the FUSE spectral 
analysis have already been published, and we will reiterate the 
results obtained from these previous studies. The FUSE obser- 
vations of 1 1 HAeBes are presented here for the first time. Some 
stars have been observed twice (e.g. AB Aur) or even three times 
(e.g. jS Pictoris, HD 104237) with a time separation between the 
FUSE observations ranging from about one month to more than 
one year. For the sake of consistency, we thus re-did the spectral 
analysis of the entire sample and took advantage of any new ob- 
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Table 1. Log of the observations. 



Star 


ff 


6 


FUSE 


Principal 


Total observing 


FUSE 


Published 




(2000) 


(2000) 


Program 


Investigator '"* 


time (ks) 


aperture 


Data (*) 


P Pictoris 


05 47 17.09 


-5103 59.45 


P219 


A. Vidal-Madjar 


24.78 


LWRS 


(1) 








Q119 


M. Deleuil 


34.50 


LWRS 










C132 


J.-C. Bouret 


47.79 


LWRS 




HD 135344 


15 15 48.44 


-37 09 16.03 


Q306 


A. Lecavelier des Etangs 


8.44 


LWRS 




HD 100453 


11 33 05.57 


-54 19 28.54 


C126 


archives 


11.79 


LWRS 




HD 36112 


05 30 27.53 


+25 19 57.08 


Q319 


M. Deleuil 


7.10 


LWRS 




HD 104237 


12 00 05.08 


-78 11 34.56 


P163 


archives 


18.97 


LWRS 










P263 


archives 


20.84 


LWRS 




HD 163296 


17 56 21.29 


-21 57 21.88 


P219 


A. Vidal-Madjar 


15.92 


LWRS 


(2) 








Q219 


A. Lecavelier des Etangs 


16.16 


LWRS 




NX PUP 


07 19 28.26 


-44 35 11.28 


Z906 


Observatory Program 


17.80 


LWRS 




AB Aur 


04 55 45.84 


+30 33 04.29 


P119 


A. Vidal-Madjar 


13.66 


LWRS 


(3) 








P219 


A. Vidal-Madjar 


15.55 


LWRS 




HD 141569 


15 49 57.75 


-03 55 16.36 


Q319 


M. Deleuil 


6.79 


LWRS 


(4) 


HD 100546 


1133 25.44 


-70 1141.24 


P119 


A. Vidal-Madjar 


10.35 


LWRS 


(5) 








P219 


A. Vidal-Madjar 


11.94 


LWRS 




HD 109573 


12 36 01.03 


-39 52 10.22 


B091 


archives 


12.38 


LWRS 




HD 176386 


19 0138.93 


-36 53 26.55 


P119 


A. Vidal-Madjar 


14.98 


LWRS 










P219 


A. Vidal-Madjar 


27.39 


LWRS 




HD 250550 


06 01 59.99 


+ 16 30 56.73 


B038 


C. Catala 


9.06 


LWRS 


(6) 


HD 85567 


09 50 28.54 


-60 58 02.96 


Z906 


Observatory Program 


7.06 


LWRS 




HD 259431 


06 33 05.19 


+ 10 19 19.98 


B038 


C. Catala 


16.68 


LWRS 


(6) 


HD 38087 


05 43 00.57 


-02 18 45.37 


A063 


archives 


31.12 


LWRS 










P116 


archives 


4.00 


MDRS 




HD 76534 


08 55 08.71 


-43 27 59.86 


B038 


C. Catala 


5.71 


LWRS 


(7) 



*°* Most of the data were obtained by our Team (through guaranteed time and guest investigator observations). Some of the data are from the 
FUSE archives or from observatory programs dedicated to instrument testing. 

'*': data already published. 

References: ( 1) Lecavelier des Etangs et al. (2001); (2) Deleuil et alj ( 120051) : (3) iRoberge et al] ( 1200 Ih : (4) iMartin-Zaidi et al.l ( 120051) : (5) 
iDeleuiletal] ( 12004) : (6) Bouret et al. (2003) : (7) Martin et al. (2004) . 



servations to improve the S/N of the spectra and the H2 spectral 
analysis. 

The fundamental parameters of the stars are given in Table |2] 
Most are from the literature, but in some cases, we re-evaluated 
the effective temperature s of the stars (for details on the 
method, see B ouret et al] l2003h . Using these values and the 
iPalla & Stahlen (|1993|) pre-main sequence tracks, we calculated 
the ages and ma sses using interpolation routines written by 
iTesti et al.1 (Il998l) . All the stars are plotted on the HR diagram 
presented in Fig.[Tl showing their evolutionary status. 

All the stars were observed with the FUSE 30" x 30" LWRS 
aperture at a resolution of about 15 000. The FUSE observations 
cover the wavelength range from 905 A to 1 187 A. For the sake 
of consistency, all the spectra were reprocessed with version 
3.0.7 of the FUSE caUbration pipeline, CalFUSE (! Dixon et alj 
l2007h . For each star, the total exposure time was split into sev- 
eral sub-exposures which have been co-aligned using a linear 
cross-correlation procedure and added segment by segment. In 
some cases, the co-added spectra in each detector channel were 
rebinned in wavelength to increase the S/N ratio without degrad- 
ing the resolution. 

3. Molecular hydrogen analysis 

3.1. H2 spectral analysis 

The FUSE spectral domain offers access to the Lyman and 
Werner series electronic transitions of H2. These electronic 



transitions are between the rotational levels (J) of the ground 
vibrational level (v - 0) of the fundamental electronic 
level (X) and the rotational levels of the vibrational lev- 
els of the first (B) and second (C) excited electronic levels. 
Under certain conditions, spectral lines coiTesponding to tran- 
sitions from the rotational levels of higher vibrational lev- 
els (v > 1) of the fundamental electronic level can also be 
observed (e.g. Lecaveli er des Etangs et al. 2003; Bour et et al] 
l2003t iBoisse etani2005h . 

For the hottest stars in our sample, typically stars whose 
Tetf values are greater than 10000 K, the FUV stellar flux is 
high enough to observe circumstellar and/or interstellar (CS/IS) 
lines in absorption against the continuum. For the coolest 
stars, the lack of sufficient flux in the continuum prevents 
easy detection of CS/IS features. However, as highlighted by 
iLecaveher des Etangs et alj (1200 lb . stellar emission lines, when 
present, can provide "continuum" flux for CS/IS absorption fea- 
tures. Most of the cool stars in our sample exhibi t broad emis- 
sion i n the Ovi /1/11032-1038 resonance doublet dPeleuil et al.l 
'2006"). A number of H2 transitions occur in the spectral domain 
near 1034 A. In particular, the strongest transitions of H2, which 
arise from the 7 = to 7 = 2 levels, fall in the same wave- 
length range as the 1038 A emission line of the Ovi doublet. 
This line is the weaker component of the O vi doublet and, as al- 
ready shown bv lRoberge et alj ( 1200 lb . the flux in the blue part of 
the emission profile can be completely suppressed by strong H2, 
but also by C 11 absorption dRoberge et al. 2006), making detec- 
tion of H2 challenging. In that case, other lines coiTesponding to 
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Table 2. Physical parameters of the sample stars. 



Name Spectral T^ff E (B-V) v sin / 
Type (K) (kms'') 



Vrad 

(kms-') 



Age (W 
(Myr) 



Mass "" 
(Mo) 



Distance 

(pc) 



LogL. 

(Lo) 



Ref. 



P Pictoris 


A5V 


8200 




0.02 


140 


+20 


20 


1.75 


19.3 


0.94 


1,2,3,4,5 


HD 135344 


F4V 


6750 




0.15 


80 


-3 


30 


1.3 


84 


0.48 


1,6,7,8,9,11 


HD 100453 


A9 


7500 




0.08 


39 


- 


10 


1.7 


114 


0.95 


1,7,8,10,11 


HD 36112 


A5 


8120 




0.07 


60 


+ 17.6 


5.0 


2.0 


200 


1.35 


1,9,11,12,13 


HD 104237 


A4 


8500 




0.10 


10 


+ 14 


2.4 


2.4 


116 


1.55 


1,7,11,12,14 


HD 163296 


AlVe 


8700 




0.06 


133 


+4 


4.7 


2.2 


122 


1.48 


1,12,15,16 


NX PUP 


AO 


9770 




0.19 


120 


- 


2.8 


2.4 


450 


1.63 


1,12,18,19,20 


AB Aur 


AOV 


9800 




0.16 


100 


+21 


2.6 


2.4 


144 


1.68 


1,23,24,25 


HD 141569 


B9V 


10040 




0.13 


236 


-6.4 


>5.0 


2.3 


108 


1.35 


1,12,28 


HD 100546 


B9V 


10470 




0.08 


55 


+ 17 


>9.0 


2.4 


103 


1.51 


1,11,12,29 


HD 109573 


A0/B9 


10500 




0.03 


152 


+9.4 


>8.0 


2.2 


67 


1.37 


1,17,21,22 


HD 176386 


B9/B8 


12000 


(a) 


0.19 


220 


+7.3 


3.5 


2.7 


140 


1.69 


1,30 


HD 250550 


B7 


12800 


(a) 


0.22 


110'"' 


+31 


1.0 


3.6 


> 160 


2.19 


1,12,23,32 


HD 85567 


B5V 


15200 


(a) 


0.23 


60 w 


0/-5 


1.0 


4.2 


>480 


2.54 


1,12,33 


HD 259431 


B5 


15900 


(a) 


0.26 


95(a) 


+43 


> 1.0 


4.4 


290-800 


2.52 


1,12,23,32 


HD 38087 


B5V 


16500 


(a) 


0.29 


100 « 


+33 


> 1.0 


4.6 


381 


2.61 


1,14 <■=' 


HD 76534 


B2 


20000 


(a) 


0.32 


110<"' 


+ 17 


>0.5 


>5 


> 160 


2.83 


1,12,23,34 



*"': Teff and v sin i determined from our m odeling (see | Bouret et al.ll2003l . for details of the method); 

*': Ages and masses computed by us from lPalla & Stahleii ( ll993h tracks using the methods described in lTesti et al.l ( I199m) : 

'■'^': Luminosity calculated using bolom etric correction. 

References: (1) SIMBAD data base; (2) Smith & Terrile ( 1984); (3) Barrado v N avascues et al. (1999); (4) Bouret et al. ( 2002); (5)'Cri fo et all 
(1997); (6 ) Malfaitetal. ( 1998 ); (7) Meeus et al. (2001); (8) Dominik et al. (2003|l]_(?l!Dent et al. (2005); (10) Meeus et al. (200^ ; (11 ) 
lAcke & W aelkens (2004); (12) Ivande n Ancker et al. ( 199J) ; (13) [g eski-ovnava et al. (1999); (14) Bertout et al. (1999); (15) Cidale etly 
(L2OOII); (16) Mora et al. (20011); (17)lTorres et al. (2003); (18)ICorcoran & Rav ( 1997); (19) Valenti et al. (2000); (20) Bohm & Catala ( 1995E 
(21) Rover et al. (2002); (22) Gerbaldi et al. ( 1999); (23) Finkenzeller & Jankovics ( 1984); (24) Bohm & Catala ( 1993); (25) Bouret & Catap 
(fl99a); (26),Lada a985.'l; (27) .Hernandez et al.. (.2004); (28),Dunkin et al. ( 1997.^; (29).Donati et aU a997.); (3q).Siebenmorgen et al.. (.2000i) : 
(3 D iMillan-Gabet et al.l mOVt ): (32) lBouret et alj ( l2003h ; (33)lMiroshnich enko et alj ( l200ll) ; (34) lMartin et alj ( l2004l) . 



transitions arising from upper J-levels should be observed, espe- 
cially, the spectral lines from transitions arising from the 7 = 3 
and 7 = 4 rotational levels, which fall in the spectral domain of 
the 1032 A O VI emission line. 

We detected and measured H2 absorption lines in the spec- 
tra of all but 6 stars. When H2 is detected, we identified ab- 
sorption lines of H2 arising from the rotational levels of the 
ground vibrational state. Transitions arising from rotational lev- 
els as high as 7 = 5 of the first excited vibrational state (v = 1) 
of the ground electronic level are also detected in some cases. 
We performed molecular and atomic absorption line analysis 
using the Owens profile fitting pro cedure written by Dr M. 
Lemoine dLemoine et al.l 12002; Hebrard et al.ll2002l) . which al- 
lowed us to simultaneously fit all lines of each species. The 
wavelengths and oscillator strengths of the H2 lines were ob- 
tained from ^ bgrall et al. (1993a) for the Lyman system and 
lAbgrall et al.l (Il993 b) for the Werner s ystem. The inverse of the 
total radiative lifetimes are reported in lAbgrall et al.l (I2000,) . For 
details on the Une analy sis procedure, we refer the reader to 
our previous studies (e.g. lBouret et al.ll2003t iMartin et al.ll2004l; 
[Martin-Zaidi et al. 2005]). We measured the column densities at 
each energy level of H2, the radial velocities, and the line widths. 
For most of the stars, the H2 mean projected velocity is very 
close to that of the star, within the resolution of the FUSE spec- 
tra (~15 kms"'). This is a clue that the detected gas is CS in ori- 
gin. Our results are given in Table|3]with 2-cr error bars. Table|4] 
available as on-line material, contains the following information. 
Columns 1 and 2 list the number of the vibrational (y) and rota- 



tional (J) level. Column 3 gives the statistical weight of each 
level. Column 4 gives the energy of each level, and Columns 5 
to 15 present the column densities of each energy level of H2, 
when detected, for all the targets stars. 

Using the column densities measured for each level, we plot- 
ted the excitation diagrams. These diagrams show the ratio of 
the column density to the statistical weight of each energy level 
versus the energy of the level. They characterize the excitation 
conditions of the gas and give information about the temperature 
of the gas and its physical conditions. For the stars in our sam- 
ple, the diagrams show marked differences between the Herbig 
Ae/B9 stars, on the one hand, and the Herbig Be stars, on the 
other. We distinguished between these two sub-groups of stars 
in our analysis (see Sects. |4| and |5]l. 

3.2. Hi Modeling 

Absorption line spectroscopy does not allow us to directly mea- 
sure the distance of the absorbing material along the line of sight. 
In a first attempt to better diagnose the origin of the detected 
molecular gas and its excitation conditions, we modeled the CS 
H2 spectra using a photodissociation regi on (PDR) model. W e 
used the PhotoDissociation Region code of lLe Petit et al.l(l2006l) . 
referred to as the "Meudon PDR Code". The code has been de- 
veloped as a tool to investigate the physical conditions in both 
diffuse and dense IS clouds. The physics included in the code, 
as well as the numerical tech niques and several re cent improve- 
ments, are fully described in lLe Petit et alj (l2006l) . We refer the 
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Table 3. Total column densities, radial velocities, intrinsic line widths b, and origin of the detected H2 towards each star in the 
sample. 



Name 



log N(H2) 



W 



P Pictoris 
HD 135344 
HD 100453 
HD 36112 
HD 104237 
HD 163296 
NX Pup 
AB Aur 
HD 141569 
HD 100546 
HD 109573 
HD 176386 
HD 250550 
HD 85567 
HD 259431 
HD 38087 
HD 76534 



V.-ad (H2) *"' 

(kms-') 



b 
(kms" 



Origin 
of Ho 



Ref. 



inclination 
of the disk O'*-* 



Ref. 



< 17.41 

< 17 

< 16 

< 17.85 
18.68!«f 

< 19.30 

20-03-fils 
20.32+«-^« 

16 46+''--' 

< 15.40 

20.80!«;'^ 

in T^+0.20 

20 64+''" 
20 43+''" 
20.64:»-^ 



1 4+0.9 
0±2 



2±3 
20 
0±2 

~-7.0 
-1 0+" 

4 5*°-^ 
13.0+f2 



5 3+" 

J.J Q 

2 2+" 



6.8+1 

5.37 
3.4+ 



5 6+" 

J.U„ 

6.0+" 



2.36+ 
0.0+^ 



4 5+' 

T.J_| 

3.7+; 

5.0+1 



- 


(1) 


~0 


(6) 


- 


(1) 


79±2 


(7) 


- 


(1) 


-65 


(8) 


~ 


(1) 


53-57 


(9) 


CS 


(1) 


78 


(10) 


CS 


(2), (1) 


60±5 


(11) 


- 


(1) 


? 




i/CS? 


(1),(4) 


55-63 


(13) 


IS 


(1) 


51+3 


(14) 


CS 


(2), (1) 


51±3 


(15) 


- 


(1) 


20.5±3 


(12) 


CS 


(1) 


n 




CS 


(3) 


n 




CS 


(1) 


n 




CS 


(3) 


n 




CS 


(1) 


n 




CS 


(5) 


n 





*"' Upper limits on the column densities have been obtained with b ■ 
"'' Radial velocities are in the stars' rest frames. 



1 kms 



*'^' The inclination angles for the disks are given with respect to the line of sight (; = 0: edge-on), "n" indicates stars with no evidence of disk. 
References in column 6 conc ern the results given in columns 2 to 5, and references in c olu mn 8 are for the inclin atio n angles of the disks . 
Re ferences : (1) this work; (2)lLecavelier des Etangs et all |l2003l); (3)lBouret et al.l(l2003l); (4)lRoberge et all |l200ll); (5)l Martin et~ZI (|2004|): 
(6)ISrnith & Terrile ( 1984); (7) Dent et al. (2005); (8) Dominik et al. (2003); (9) Eisner et al. (2004); (10) Gradv et al. (2004); (ID lGradv et al.l 
(I2000h ; (12) iAugereau et al.i ( il999.) ; (13) iEisner et alj ( i2003.) ; (14) iWeinberger et al.,(il999i) ; (15) iAugereau et al.i ( i2001.) . 



reader to this paper and provide only the main outlines of the 
code here. 

The code models a slab of IS medium, in a stationary state, 
heated by one star's radiation field and the galactic field in a 
plane-parallel geometry. The IS radiation field flux is assumed 
isotropic on both sides of the cloud, while the added incident 
stellar radiation is normal on on e side. The code uses a standard 
galactic field law (lDrainel[l978l) . The stellar FUV radiation field 
is given by a blackbody law calculated from the star's radius 
and temperature and the star-to-cloud distance. The model input 
parameters are thus the radius and effective temperature of the 
star, the distance of the observed gas from the star, the density of 
the medium, and the extinction (Ay). The equations of thermal 
and chemical balance as a function of depth within the cloud are 
solved in a plane -parallel geometry with 113 chemical species 
and 967 chemical reactions included. As a result, the code pro- 
vides the temperature distribution within the gaseous medium 
and the abundances of the chemical species derived from the 
chemical balance. 

Without any precise information on the nature of the dust 
grains in the absorbing medium, we used s tandard ISM value s 
for the grain size distribution and properties (iMathis et al.l 19771) . 
When available, we used lUE spectra to determine the column 
density of H i, which allowed us to estimate the ratio of the gas 
column to the Ay as an input parameter. We used the standard 
ISM value of this ratio when no observational constraints were 
available. For most of our targets, little or no information is avail- 
able in the literature about the properties of the observed gas, 
namely its distance to the star or its density and temperature. We 
thus proceeded by an iterative method, varying both the distance 
and density of the gas to obtain the best fits to the observed ex- 
citation diagrams. We checked that the model output parameters 



are consistent with the molecular fraction we derived from our 
spectral analysis (see Sect. |5]l. We also compared the synthetic 
spectra provided by the code to the observed ones, as an addi- 
tional check. The results of our modeling are presented in the 
next sections. 



4. Ae/BS stars 

In this section, we first present the stars for which we have not 
detected H2 absorption lines. In the following subsections, when 
H2 is detected, we distinguished the stars for which we observed 
cold H2 and warm H2 on the basis of their excitation diagrams. 

4.t. Stars lacking H2 absorption iines 

We could not detect any H2 absorption lines for 6 stars in our 
sample. 

jS Pictoris: Early FUSE observations of/? Pictoris in 2000 
allowed us to report a clear deficiency of molecular hydrogen in 
the d isk of the star, consistent wit h the evolved status of the sys- 
tem dLecavelier des Etangs et al]|2001 ). Taking advantage of a 
third exposure acquired in November 2002 with a total exposure 
time of about 50 000 seconds, much longer than the two previ- 
ous ones, we checked for any H2 absorption features in the new 
spectra (Fig.|2]i. Co-adding all three spectra allowed us to refine 
the upper limit on the total H2 column density previously esti - 
mated at about 10^^ cm"^ bv lLecavelier des Etangs et al](l2001h . 
Assuming a low value for the intrinsic line width (b = I km s~'), 
we found the new estimate to be less than 2.6x10'^ cm"^. 

HD 109573: HD 109573 is the only star in our study for 
which no H2 absorption lines could be detected, despite a notice- 
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Fig. 2. Ovi emission doublet near 1032 A and 1038 A in the 
FUSE spectrum of/? Pictoris observed in 2002. Here, only the 
exposures obtained during the night have been co-added, to 
avoid contamination by strong airglow emission lines due to 
the Earth's atmosphere. The positions of the H2 transitions we 
should observe if H2 was present along the line of sight are plot- 
ted. The H2 levels are mai-ked H20 for (v = 0,7 = 0), H21 
for (v = ^^J - 1), etc. For details on the Cii absorption, see 
iRobergeet all (120061) . 



able stellar continuum flux. The stellar flux is well-developed in 
the reddest part of the FUSE wavelength range, above 1 1 10 A. 
Below this wavelength, where the strongest H2 transitions fall, 
the stellar flux level is much lower but would be high enough 
to observe H2 absorption lines if present, but the signal-to-noise 
ratio is very low (S/N ^ 1.5) and no stellar emission lines are 
present to overcome this issue. 

HD 109573 is also an evolved system, the star being very 
close to the main sequence. It is known to possess a flat debris 
disk seen nearly edge-o n, with an inclination a ngle of about 20° 
from the line of sight (lAugereau et al.|[l999l) . The dusty disk 
has two components: an inner annulus that extends between 
4-9 AU to ~55 AU from the central star, and an outer ring at 
70 AU (Wahhaj et al. 2005; Augereau et al. 1999). These two 
components are separated by a gap located between 55 and 60 
AU from t he star, which is proba bly the consequence of planet 
formation dAugereau et al.l [l999l) . The outer ring is about 15 
AU th ick, with an inclination angle of 13±1° from the line of 
sight JTelesco et al . 2000). This implies that the upper layers 
of the disk intersect the line of sight and could potentially 
produce a bsorption features. However, as previously reported by 
IChen & Kamp (2004), we do not observe absorption lines of H2 
in the FUSE spectrum. We set an upper limit on the total column 
density of H2 of < 2.5 x 10'^ ^"^^^ w hich is consistent with 
the limit found bv lChen & KampI (|2004|) . This non-detection of 
H2 implies that the disk of HD 109573 is almost depleted of 
molecular gas, which agrees with the evolutionary status of this 
star. 

NX Pup, HD 36112, HD 135344, and HD 100453: The 

photospheric continuum of these four stars is too low in the FUV 
(S/N ~2 at 1 100 A) for a reliable detection of CS/IS absorption 
lines. However, as in the case of /J Pictoris, we took advantage 
of the O VI emission lines to place limits on the total H2 col- 



umn densities along the line of sight. The presence of a CS disk 
around NX Pup is suspected from its SED. For the three other 
stars, all younger than /? Pictoris, the non-detection of H2 is con- 
sistent with the high inclination angles for their disks, which do 
not intersect the lines of sight to the stars (see Table|2]l. It shows 
that, if present, the CS gas does not extend far enough above the 
midplanes of the disks to be detected in absorption and that the 
remnant of the initial CS envelopes have been dispersed or are 
too low in density to allow detection. 

4.2. ColdH2 

HD 141569: The HD 141569 disk is usually called a transitional 
disk, in which there are signs of dust clearing in the inner 
disk, but significant primordial molecular gas is still supposed 
to be present in the disk . In a previous detailed analysis (see 
iMartin-Zai'di et al.] l2005h . we demonstrated that the absorp- 
tion lines of molecular hydrogen and several atomic species 
observed in its FUSE spectrum are very likely related to the 
slightly reddened dif fuse outer region s urrounding the L134N 
dark cloud complex (iJuvela et al.ll2002h . The non-detection of 
CS gas implies that there is no remnant CS envelope and that 
all the gas has had time to collapse into a flat or very slightly 
flared disk. This result is co nsistent with the ^^ CO observations 
at 345.796 GHz described in lDent et al.1 (l2005h . 

AB Aurigae: The presence of H2 in the far- UV spectra of 
ABAurigae was first reported and analyzed by iRoberge et al.l 
(1200 Ih from an early FUSE observation. A new observation ob- 
tained more than one year later allows us to confirm and refine 
the analysis. 

We identified and measured H2 lines arising from 7 = to 
7 = 6 in the ground vibrational state (v - 0). As suggested by 
iRoberge e t al. (2001), the H2 could be IS in origin or have a CS 
origin, as a cold rem nant of its natal cloud or a CS envelope 
dSemenov et al.ll2005l) . The relatively low spectral resolution of 
FUSE prevents us from drawing firm conclusions about the ori- 
gin of the detected gas on the basis of its radial velocity and line 
width. The excitation diagram derived from the column densities 
in the J-levels may provide additional clues for further constrain- 
ing the location of the observed gas. For AB Aur, it confirms that 
the H2 is thermalized up to 7 = 2 with a low kinetic tempera- 
ture of 56±4 K (Fig. [3]), results in good agreement with those of 
Roberge et al. (2001). Since the relative populations of the lower 
H2 energy levels are determined primarily by thermal collisions, 
the excitation temperature found from these levels should be 
close to the kinetic temperature of the gas. On the other hand, 
the column densities of the higher J-levels are compatible with 
thermal equilibrium at a temperature of 373+83 K. 

Using the Meudon PDR Code, we checked that such exci- 
tation conditions could not be reproduced by a single absorbing 
component along the line of sight towards the star. We varied the 
different parameters within their error bars and tried to model the 
observations with different Ay values found in the literature, but 
no PDR model fit the data. There are still some free parameters 
to explore, but we need additional observational constaints to 
make an accurate modeling. As shown in Fig.|3] a diffuse com- 
ponent (n//~250 cm""*) located at 0.3 pc from the star reproduces 
the column densities of the three first energy levels, but failed to 
reproduce the higher J-level populations. 

Two-temperature behavior in excitation diagrams for H2 
is o ften observed in diffuse an d translu cent IS clouds (e.g . 
Spi tzer & Cochranll973l:rGrv et a l. 2002; B rowning et al.l2003h . 
but its origin is not understood well and remains a subject of con- 
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Fig. 3. Excitation diagram for H2 towards AB Aurigae. The ob- 
served level column densities (filled circles) show that H2 is ther- 
malized up to 7 = 2 with a low kinetic temperature of about 
56 K, while the column densities of the higher J-levels are con- 
sistent with a temperature of about 373 K. The populations of 
the first three energy levels are reproduced well by a PDR model 
(open square) with a low density of about 250 cm""* and a dis- 
tance between the star and the gas of about 0.3 pc, with an exci- 
tation temperature of about 55 K. Our model does not reproduce 
the higher energy levels' excitation. This may be due to the pres- 
ence of a second warm/hot gaseous component along our line of 
sight (see text). 



troversy, as reviewed and discussed in ISnow & McCall (l2006h . 
Several excitation mechanisms could be invok ed separa tely or 
all together to explain this behavior. IShuU & Beckwith ( Il982h 
initially suggested that the populations of the higher levels can 
be inflated by UV and/or formation pumping and radiative cas- 
cades. We checked that, to populate the J > 2 levels by UV 
pumping alone, we would have to increase the radiation field 
(G~1.5 in Draine's units) by a factor of more than 10. An alter- 
native model invokes local heating of the gas within the cloud. 
Such heating could be generated by magneto-hydrodynamic 
shocks (see Chieze et al. 1998) or transient small-scale turbu- 
lence that results in localized tiny warm regions, as suggested by 
iFalgarone & Pugetl (Il995h . 

The simplest explanation is, however, the presence of dis- 
tinct absorbing media along the Une of sight with very differ- 
ent excitation conditions (Browning et al.l2003h . For all the stars 
in this study, we found i>- values much higher than expected for 
pure thermal broadening. This is an unusual situation in diffuse 
ISM clouds, where ty pical /7-values are about 1.5 up to 3 kms ' 
(e.g. iGrv et al]|2002h. However, higher values are also rou tinely 



reported (e.g. lLehner et alll2003l; iRedfield & Linskvll2004 . The 
simplest explanation for this phenomenon is the presence of sev- 
eral gaseous components along the line of sight that could appear 
blended in the spectra. But it could also be related to gas motions 
on different spatial scales as suggested by different studies. Such 
motions could originate in turbulent and wa rm layers within the 
cloud, through tran sient shocks or vortices (iJoulain et al.lll998t 
iGredel et al.ll2002|) . which lead to increased observed /j-values. 
Without observations at higher spectral resolution or additional 
constraints on the chemistry of the gaseous medium, we cannot 
distinguish between explanations. 



In conclusion, we cannot distinguish with certainty yet 
whether the H2 gas in this source has a circumstellar or inter- 
stellar origin, but the cold component appears to be located at a 
relatively large distance from the star (~0.3 pc). 

Now the question arises, whether the gas has a circumstel- 
lar origin, i.e. could it be part of the disk or must it be part of 
a separate, more spherically distributed, components of CS gas? 
AB Aur is known to harbor an extended CS disk with an in- 
clination angle between 27° and 35° from the plane of the sky 
(lEisner et al.ll2003l: iPantin et al.ll2005h . i.e. nearly face-on. Since 
a detection of H2 in absorption means that the line of sight to- 
ward the star must pass through the H2 gas, the near-face-on 
inclination appears to make it unlikely that the gas is part of the 
disk. However, a more quantitative analysis is required to make 
any firm statements. Suppose the disk is at its thickest (geomet- 
rically) at a radius r, where it has a gas surface density Sgas and 
a gas temperature T. Assuming vertical hydrostatic equilibrium, 
the gas density as a function of height z above the midplane at 
that radius is 






(1) 



p/ 



with 



H„ 



kTr^ 



( fiMpGMt 



(2) 



where fi - 2.3 is the mean molecular weight of the molecular- 
hydrogen-helium mixture, nip the proton mass, G the gravita- 
tional constant, and k the Boltzmann constant. The column depth 
along the line of sight passing through (r, z), i.e. with the disk at 
an inclination of / = arctan(r/z) away from face-on, is then esti- 
mated as 



A^H, = 



Pgas' 



Pgas"p 



linip ^inip 



(3) 



Here we simply took as radial extent / the pressure scale height 
I ^ Hp.ln reality this may be a factor of a few higher 

Let us tak e r to be the outer rad ius of the disk estimated at 
'"out =800 AU JSemenov et al.ll2005h and a disk mass of at most 
1 Mq. With Mdisk - 2gas('"out)7r?"^m, this gives a surface density 
of about Sgas ^ 4 g/cm^ at r = rout. Using these numbers and 
the stellar parameters from Table|2] we find Hp{roudlfoM - 0.36 
for a temperature of T = 95 K, which, with an inclination of 35° 
away from face-on, yields a column density of log A^(H2) = 20 
as given in Table |3] In other words: we need a temperature of 
about 95 K to puff the disk up enough for sufficient matter to 
cross the line of sight at around 800 AU distance. Considering 
the roughness of this estimate and the closeness to the excitation 
temperature of 56 K for the low-lying levels, it shows that it is 
possible in principle that this gas is part of the outer part of a 
massive disk. But the same derivation yields a volume density of 
hh, - 3.4 X 10'*, which is much higher than the ~ 250 derived 
from the PDR model. Instead of 800 AU, we could also use 0.3 
pc s 61900 AU as reference radius r. Then any estimate would 
yield the disk as gravitationally unbound, and we cannot speak 
of a "disk". Thus, if the cool material is indeed located at such 
a distance, it may still be associated with the star AB Aurigae, 
but it would be meaningless to speak of disk material there and 
is likely to be either a remnant envelope or outflowing material 
from the disk. 
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4.3. Warm H2 gas 

We re-did the spectral ar ialyses for HP 100546 and HP 163 296, 
which were presented in lLecavelier de s Etangs et al] ( l2003l) . and 
found the same results. We also refined the analysis of the F USE 
spectrum of HP 104237 presented in He rczeg et alj (l2003h and 
found compatible results. In the HP 100546 spectrum, lines aris- 
ing from the J-levels of the ground vibrational level (v - 0) up 
to 7 = 10 and J-levels of the first vibrational level (v = 1) up to 
/ = 5 are detected. In the HP 163296 and HP 104237 spectra, 
only lines arising from the J-levels up to v = 0, 7 = 4 and v = 0, 
J - 5, respectively, are observed. Higher rotational J-levels and 
vibrational levels are beyond the detection limit, because of the 
sharp flux decrease with decreasing wavelength caused by the 
later spectral types of these two stars and because of the lower 
signal-to-noise ratios of the data. These low signal-to-noise ra- 
tios are responsible for the rather large uncertainties in the deter- 
minations of the column densities (Fig.|4]i . 

We plotted the excitation diagrams of H2 for each of the 
three stars, and found two-temperature behavior in the case of 
HP 100546 (Fig.EJi. In contrast to the sources in Sect.|4|2] for 
all three stars, the H2 is thermalized up to 7 = 4 with surprisingly 
high kinetic temperatures, greater than 300 K and as high as 
758 K in the case of HP 100546 (Fig.|4]i. In addition, the temper- 
ature derived for the higher J-levels observed in HP 100546, is 
1495±70 K. Such excitation conditions clearly differ from those 
observed in the interstellar medium. In addition, the measured 
radial velocities favor a CS origin for the H2 at the FUSE spec- 
tral resolution. 

Us ing the numerical model developed by iLe Bourlot et al.l 
( Il993h . we estimated that the warm H2 (with 7 < 4) should 
be located at a distance of about 1.5 AU from HP 100546, 
4 AU from HP 163296, and 9.5 AU from HP 104237 to ex- 
plain the obs erved excitation conditions (for d etails about the 
method see iLecaveher des Etangs et al.l 12003"). The observed 
temperature cannot be explained by the Galactic UV radiation 
field, which implies that the observed H2 is circumstellar and 
lies very close to its central star. Using the volume densities de - 
rived from the critical density given bv lLe Bourlot et all (Il999h . 
we found typical thicknesses of <0.2 AU for the absorbing layer 
towards HP 100546, <1.2 AU for HP 163296 and 3.2 AU for 
HP 104237. These results reinforce our interpretation of a CS 
origin for the detected gas. 

This again poses the question of whether this material is 
part of the disk or part of some more spherically distributed 
gaseous circumstellar material. For all three sources, the incli- 
nation is more than 50° away from edge-on. Using Eq. [3] with 
a gas surface density of, say, Sgas = 10"^ g/cm^, at the radii es- 
timated above, and inclination of 50° and required column den- 
sities given in Table |2l we find that we need a temperature of 
1 1600 K for HP 100546, 4700 K for HP 163296, and 2300 K for 
HP 104237 to support hydrostatic equilibrium and yet have the 
required column density of H2. All these temperatures are well 
above the excitation temperatures measured for these sources, 
so the disk is excluded as the carrier of the warm H2 gas we 
detected. 

We can think of two possibilities to explain the excitation 
of this gas. The first is that this medium is related to a non- 
radiatively heated medium close to the stellar surface, such as 
a chromosphere, as prove d by the well-known FUV and X-ray 
activity of these stars (e.g. 'Peleu il et alj2004l2005l:lGradv et"al] 
[2004). Others Herbig Ae stars exh ibit noticeable activ ity in the 
far-UV but at a much lower level dPeleuil et alj|2006l) . The H2 
we observe could be produced in the very outer layer of this ex- 
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Fig. 4. (a) Excitation diagram of H2 for HP 100546. The H2 
is thermalized up to / = 4 with a kinetic temperature of 
758+147 K. The column densities of the higher J-levels are con- 
sistent with a temperature of about 1495+70 K. (b) Excitation 
diagram of H2 for HP 163296. As for HP 100546, the H2 is 
thermalized up to 7 = 4 with a lower kinetic temperature of 
432+135 K. Here, only the v = levels up to 7 = 4 are observed 
because of the later type of the star (Al Ve) and the low S/N ratio 
of the spectrum, (c) Excitation diagram of H2 for HP 104237. 
The H2 is thermalized up to 7 = 4 with a kinetic temperature of 
306±80 K. The higher J-levels may be consistent with a higher 
temperature, as suggested by the 7 = 5 column density. 



tended hot region. Pue to the inclination angle of the disk, this 
region should have a more or less spherical geometry to allow 
gas to be observed in absorption. 
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Another natural origin for the excited H2 is photoevapo- 
ration of the disk by FUV stellar radiation (Hollenbach pri- 
vate communication). This is a mechanism of disk dissipation 
similar to the better kn own EUV-driven photoevaporation (e.g. 
iHollenbach et al.l Il994t) . but this time driven by non-ionizing 
FUV photons. In contrast to the ionized EUV-driven wind, which 
originates as close as ~1 AU from the star, the FUV-driven wind 
is neutral and originates in the outer regions of the disk (many 
tens of AU away from the star). It is expected to have a much 
stronger wind mass loss rate than the EUV wind. According to 
Gorti & Hollenbach (ApJ submitted), such strong disk winds are 
expected to flow away from the source in a nearly spherical way, 
and they are therefore interesting candidates for explaining the 
observed warm circumstellar matter The problem in this case is, 
however, that this FUV driven wind originates much farther out 
than the few AU radii derived from the excitation diagrams. 

Further observations, such as IR spectra of H2, are required 
to give a better understanding of the origin of the gas and the 
physical processes that produce this excited H2. More modeling 
of photoevaporating gas or chromosphere models, which would 
take non-LTE excitation mechanisms into account, could also 
help to better constrain the excitation conditions. 



Table 5. Column densities of H i measured from the lUE spectra, 
when good data were available. 



Stars 



Column density 
of H I (cm"-) 



N(H,„,)/E(B-V) 
(atoms cnT^ mag"') 



/ 



HD 250550 
HD 259431 
HD 38087 
HD 76534 
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of the excitation temperature and molecular fraction. The best- 
fitting model allowed us to derive the total density and the dis- 
tance of the gas from the central star, which are not constrained 
by absorption lines observations. Our results are tabulated in 
Tablei] 

Table 6. Distance of the gas from the central star and total den- 
sity of the observed medium, as determined by the best-fitting 
model using the Meudon PDR Code. 



5. Be stars 

This second group of stars contains all the stars earlier than B9 
type, which are also the youngest stars in the sample. For all but 
one, namely HD 250550, we have identified and measured lines 
arising from the rotational levels of the ground vibrational level, 
up to at least J - 1, and also lines from the first vibrational 
level. The derived excitation diagrams, plotted in Fig. |5] show 
that the H2 is thermalized up to 7 = 3 with a kinetic temperature 
around 100 K. The temperatures given by the column densities 
of the high J-levels range from ~500 K to ~1600 K. The radial 
velocities of this matter suggests that it is of circumstellar origin. 
As seen in Fig.|5] the excitation diagrams point to similar ex- 
citation conditions for all the stars, which is a clue to a common 
structure in their CS environments. Moreover, since the v sin / of 
these stars span a wide range of values (see Table |2| and since 
the detection of H2 requires the line-of-sight to go through the 
material, it is unlikely that the material in all these objects is part 
of a circumstellar disk. Instead it is more likely that this matter 
is part of the remnants of the clouds within which the stars were 
formed. 



5.1. Modeling the H2 excitation 

Little is known about the CS environments of these young stars 
and to model the physical conditions of the gas we detected, we 
had to estimate realistic values for the input parameters. To de- 
rive the H I column density for each star, we used archival lUE 
spectra to model the Lya line (Table |5]l, when good data were 
available. Using the E(B-V) values listed in Table |2] we calcu- 
lated the gas-to-dust ratios N(H,o/)/E(B-V), in the environments 
of the stars and the molecular fractions (/), i.e. the fraction of 
hydrogen atoms in molecular form (Table |5]). Our gas-to-dust 
ratios values are comparable to the mean interstellar value of 
about 5.8 X 10^' atoms cm"^ ma g"\ as measured fo r standard IS 
clouds within 2 kpc of the Sun dBohlin et al.lll978h . This favors 
an IS origin for the H2 we observed towards these stars. 

For each star, we run grids of models using the Meudon PDR 
Code, with different distances and densities for the gas, to find 
the best-fitting model that is consistent with the observed values 



Star 


Gas to star 


Gas volume 




distance (pc) 


density, Hh (cm"^^) 


HD 176386 


0.04 


3000 


HD 250550 


0.03 


600 


HD 85567 


0.10 


250 


HD 259431 


0.35 


70 


HD 38087 


1.60 


12 


HD 76534 


0.40 


170 



As shown in Fig. |6l the general shape of the excitation dia- 
grams is reproduced well by PDR models with a single compo- 
nent, but some details remain. The temperatures on the plots cor- 
respond to the excitation temperatures of the low energy levels 
(7 = 0-3) given by the best-fitting models. These temperatures 
are fully consistent with those observed, and always fall within 
the error bars of the observed excitation temperatures (Fig. |5]). A 
relatively diffuse and cold component reproduces the excitation 
of the lower J-levels, which includes more than 90% of the de- 
tected gas. However, we emphasize that, as for the HAe stars, our 
modeling fails to reproduce perfectly the excitation of the first 
J-levels of the ground vibrational level and the higher J-levels 
(7 > 6, V = and J, v - 1) simultaneously. By increasing the 
FUV field, one could expect to increase the excitation of these 
levels, but this would affect the column densities of the lower 
J-levels due to photodissociation. This shows that UV pumping 
and formation pumping are probably not the only mechanisms 
playing a role in the excitation of intermediate J-levels we ob- 
serve (see Sect. [3]). However, at the present time, the H2 pumping 
formation process is not completely understood, and one could 
expect that a slightly different modeling could vary the fit. 

Using the photospheric models (see Sect.|2]and iBouret et al.l 
2003) and the attenuation of the UV continuum provided by the 
PDR code, we computed the synthetic spectrum of each star. To 
calculate the spectrum, we only took the first three J-levels of 
H2 into account, which produce the strongest absorption lines. 
Figure I2] shows the synthetic spectrum of HD 76534 overplot- 
ted on the observed FUSE spectrum. The synthetic spectra of 
the other stars are available on-line (Fig. [8] to Fig. [T2] |. The 
good agreement between the modeled and the observed spectra 
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Fig. 5. Excitation diagram for H2 towards the Be stars in our sample. For each star, the H2 is thermalized up to / > 3 with a kinetic 
temperature of about ~ 100 K. 
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Fig. 6. Best-fitting models obtained with the Meudon PDR Code overplotted on the excitation diagrams derived from the observed 
column densities (Fig.|5]l. The general shape of the diagrams are well-reproduced, but the excitation conditions of the higher J-levels 
cannot be reproduced well by the model (see text). 
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Fig. 7. Overplot of the FUSE spectrum of HD 76534 (thin line) and the synthetic spectrum obtained by combining our photospheric 
model with the results of the Meudon PDR Code (thick line). This plot only presents the portions of the spectrum where the stronger 
absorption lines of H2 (v = 0, 7 = 0, 1 , 2) are observed. These principal prominent absorption troughs caused by molecular hydrogen 
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strengthens confidence in our estimates of the distance, density, 
molecular fraction, and gas-to-dust ratio we estimated for each 
star from the different observational data. 



5.2. Additional information on tine targets 

- HD 176386 - This star i s the primary s tar of a binary 
system (Proust et all 119811: IVaz et all I199&) . Our analysis 
of the FUV spectrum of the star is consistent with the star 
being surrounded by a large CS envelope (d = 0.04 pc), 
a result in agreement with previous studies in different 
spectral domains. ISiebenmorgen et al.l (I2OOOI) showed that 
the FWHM derived from the ISOCAM observations are 
indicative of a large extended halo around HD 176386, 
which was confirmed by the ISOPHOT multi-aperture 
sequence at 7.3/jm. 



located in the nebula NGC 2024, but the resolution of our 
observations does not allow us to rule out the presence of a 
second component corresponding to the star's envelope/disk. 



- HD 250550 - In iBouret et al.1 (l2003b . we showed that the 
properties of the H2 and the atomic species derived from 
the FUSE spectrum of HD 250550 indicate that we are 
probing a dense CS environment, related to the remnant of 
the molecular cloud that collapsed to form the star. We refer 
the reader to this paper for more details. 

- HD 85567 - Although our analysis does not allow us to rule 
out the presence of a CS disk, it does favor an interpretation 
in terms of a large CS envelope (li = 0.1 pc). The presence 
of two components with different temperat ures in the 
excitation diagram agrees with the Malfait e t al.l (Il998l) 
study, which fit the SED of this star with a two-component, 
optically-thin, dusty envelope model. The apparent dip in 
the SED of the star between 6 and 1 fim was interpreted as a 
physical hole in the dust distribution, caused by the break-up 
of an optically thin dusty disk (Lada & Adams 1992]). In 
addition, iMiroshnichenko et al.l (12001 ) have shown that, 
without imaging observations, it is impossible to distinguish 
between an extended, optically-thin, mostly spherical 
component, and a more compact disk. Those authors show 
that the emission-line profiles indicate complex structure 
in the CS envelope, which is likely non-spherical, with an 
optically-thick component and the possible presence of an 
extended optically-thin dusty component. We emphasize 
that if two components are present, we cannot distinguish 
them at the FUSE resolution. 

- HD 259431 - Our analysis, presented in a previous paper 
dBouret et al.ll2003l) . indicates that the material we observe is 
most likely related to a large CS envelop e. This is consistent 
with the study by iMalfait et all ( Il998l) . who fit the SED 
of this star with a two-component, optically-thin, dusty 
envelope model. Although the mid-IR observations of 
HD 259431 by Polomski et al. (2002) are consistent with 
a model of a moderately flared CS disk, there is no clear 
observational evidence of a disk around this star. 



HD 38087 - ISnow & Witj ( Il989l) concluded from their 
analysis of the lUE spectrum that the star is surrounded by a 
dense envelope in which unusual grain growth h as occ urred. 
These results were confirmed by iBurgh et al] (12002 ). who 
show that the region surrounding HD 38087 is dense and that 
the dust grains may be larger than average in the large nebula 
(NGC 2024) in which the star is located. From our model- 
ing, we find a very low gas density and a large gas-to-star 
distance. The H2 we observe in the FUSE spectrum is likely 



HD 76534 - From a previous analysis dMartin et al.ll2004l) . 
we concluded that, in the hostile CS environment of 
HD 76534, which is a B2 star, the presenc e of a CS disk is 
very unlikely. This agrees with the Hillenbr and et al.l (Il992l) 
analysis of the SED of HD 76534, who stressed the sim- 
ilarity of the very low NIR excess with that of a classi- 
cal Be star, which is generally thought to be due to free- 
free emission in an ionize d envelope rather than to CS dust 
dHamann & Perssonll 19921) . 



6. Evolutionary trends 

For each group of stars, we explored whether evolutionary trends 
could be found in our results, especially between the total col- 
umn densities of H2 in the circumstellar environment of the stars 
and the ages of the stars. 



6.1. Ae/B9 Stars 

As we discussed in Sect.|4] the HAe stars form a very inhomoge- 
neous group. We plotted the column densities of H2 towards the 
HAe stars versus the inclination angles of their disks in Fig. [T3] 
and in Fig. [14] the column densities of H2 versus the ages of 
the HAes of our sample are presented. Both Figs. [T3] and [T4l are 
consistent with CS material that is not in disks. They are neither 
consistent with non-associated interstellar clouds nor explained 
by absorption in the disk material. 



T 



HDl 00546 



T 



20 40 60 80 

inclinotion ongle for the disk from the line of sight (degrees) 

Fig. 13. Total H2 column densities as a function of the inclina- 
tion angles away from edge-on of the CS disks given in Table|3] 



On the other hand, when excluding j3 Pictoris and 
HD 135344, we find that the amount of H2 is higher for the 
younger stars than for the older ones (Fig. \T4\ . Such a correla- 
tion is very unexpected if we observe non-CS media. However, 
our results do not allow us to rule out CS envelopes around these 
stars. Indeed, for some of these stars, the pre sence of both a disk 
and an envelope has been demonstrated (e.g. lGradv et al.|]2001h . 
This would confirm the canonical view of star formation, which 
predicts that the envelope becomes progressively less important 
for obscuring the star and for dominating the system mass with 
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Fig. 14. Total H2 column densities as a function of the ages of 
those stars known to possess CS disks. Except for /3 Pictoris 
and HD 135344, for which we cannot conclude because we 
only have upper limits, the amount of H2 is greater towards the 
youngest stars. The correlation between the age and the amount 
of gas is shown by the dashed line. 

increasing age. The envelope dominates the mass of the system 
in the early stages of stellar formation, but when the star is close 
to the ZAMS, the envelope is expected to have dissipated be- 
cause of the action of s tellar winds prior to central clearing o f 
the disk (see reviews bv lCalvet et"aLl(2000l; iMundv et al.ll2000l) . 
The correlation between the ages of the stars and the total 
amount of H2 needs to be confirmed. For most stars, we only 
estimated upper limits on the column densities of H2, and there 
are large uncertainties on the ages of some of them . Only precise 
measurements of the amount of H2 and the ages of the stars could 
help to confirm or deny the suspected correlation. 

6.2. Be stars 

Our analysis of the FUSE spectra of the Herbig Be stars shows 
that this set of stars presents similar H2 physical properties, fa- 
voring an interpretation in terms of more or less spherical me- 
dia around these stars, very likely related to the original clouds 
in which the stars formed. We find, however, no correlation be- 
tween the total column densities of H2 and the age of the stars 
(Fig.fTSll. as one would expect from rapid dissipation of the orig- 
inal envelope. This lack of correlation could be due to different 
phenomena: (i) large contamination by the surrounding nebula 
and thus the presence of different unresolved components along 
the line of sight or (ii) uncertainty about the ages of the stars. 
Herbig Be stars are more massive and more luminous than Ae 
stars, and the time-scale to reach the ZAMS is very short, mak- 
ing their ages very difficult to estimate. 

7. Summary and conclusions 

We have distinguished two groups of stars: the stars known to 
possess CS disks (Ae/B9 stars) and the stars with no evidence of 
disks (Be stars). We interpreted the different excitation charac- 
teristics of H2 as direct probes of the origin of the observed gas: 
IS, CS disks, or CS envelopes. 

For most of the Herbig Ae/B9 stars, the incUnation angles 
of the CS disks from the lines of sight are quite high and not 
favorable to probe the disks using absorption spectroscopy. Beta 
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Fig. 15. The total H2 column densities as a function of the ages 
of the Be stars of our sample, which present no evidence of CS 
disks. 



Pictoris and the transitional object HD 109573 are exceptions. In 
the case of yS Pictoris w e confirm the non-detection of H2 in its 
FUSE spectrum given in lLecavelier des Etangs et al] (|2001|) . We 
also did not detect H2 in the spectrum of HD 109573. These non- 
detections are consistent with the evolved status of these systems 
and clearing of CS molecular gas from their disks. 

AB Aurigae is clearly not edge-on, while H2 is clearly de- 
tected, but the cool component of the gas must be located at 
such a large distance from the star that, even with a relatively 
large inclination away from edge-on, some disk gas in the very 
outer regions of the disk may still be in the line of sight. 

In the sources HD 36112, HD 135344, and HD 100453, no 
H2 is observed. This is consistent with the high inclination an- 
gles of their disks. The non-detections of gas imply that there are 
no detectable remnant CS envelopes. 

When molecular hydrogen gas is observed in the FUSE spec- 
tra of the Herbig Ae/B9 stars of our sample, our analysis shows 
several kinds of excitation, implying different origins and differ- 
ent physical processes. 

In one case, namely HD 141569, the origin of the detected 
gas is clearly interstellar. For AB Aur there is also some evidence 
of the interstellar origin on the gas, but we lack the spectral reso- 
lution that could allow us to probe the velocity distribution along 
the line of sight and definitively decide between a CS or IS ori- 
gin for the gas. In spite of its relatively high inclination away 
from edge-on, it is not geometrically impossible for the gas to 
originate in the upper layers of the outer disk regions at about 
800 AU, but the excitation conditions of H2 towards AB Aur 
are fairly typical of what is generally observed in the diffuse in- 
terstellar medium. However, we showed that the excitation di- 
agram could not be reproduced well by a single gaseous com- 
ponent along the line of sight. The large line width value favors 
the presence of two components along the line of sight, which 
are not resolved in our spectrum. In that case, one of these com- 
ponents could correspond to the remnant of the molecular clou d 
in which the star formed, as suggested bv lRoberge et al.l ( 1200 lb . 
The other component should be hotter and probably closer to the 
star to explain the excitation of the high J-levels of H2. 

For HD 100546, HD 163296, and HD 104237, we observed 
excited and probably warm/hot, circumstellar H2 that has ex- 
citation conditions clearly different from those observed in the 
diffuse interstellar medium. Such excitation conditions for H2 
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give evidence of coUisionally excited media close to the stars. 
The high values of the b parameters support this interpretation. 
In addition, the measured radial velocities favor a CS origin for 
the H2. However, assuming that the gas and dust are coupled, the 
lines of sight towards these three stars do not pass through their 
disks, and thus the H2 we observed is not located in the disks. 
This raises questions about the origin of the detected gas. One 
interesting possibility could be that the H2 gas is an FUV-driven 
photoevaporative wind from the outer parts of the disk, but we 
are still missing clues needed for a firm conclusion. The Meudon 
PDR code did not allow us to reproduce the excitation of the ob- 
served H2 gas. This implies that peculiar excitation processes 
such as shocks or X-Rays, not taken into account in the Meudon 
PDR Code, may play a role in the environment of these stars. 

The excitation conditions of circumstellar H2 around the 
stars of our second group (Be stars) are clearly different from 
those of the first group. We found similar excitation conditions 
for the H2 from one star to the next. They present similarities 
to the conditions found in the IS medium. Our analysis favors 
an interpretation in terms of spherically symmetric media, not 
affected by inclination effects. The excitation diagrams are re- 
produced nicely by PDR models, at least for the cold component 
that includes more than 90% of the gas. On the other hand, these 
models do not simultaneously reproduce the excitation condi- 
tions of the higher J-levels of H2. We very likely observe com- 
plex environments, such as large CS envelopes, remnants of the 
original clouds in which the stars formed. However, other exci- 
tation processes, in addition to those taken into account in the 
Meudon PDR Code, are probably necessary to fully explain the 
observed excitation diagrams. 

Our FUV spectral analysis reinforces the differences be- 
tween the two subclasses of st ars (Ae and Be) already high- 
lighted by different authors (e.g. lNatta et al.ll2()00l) . It also points 
to the need for complementary observations to better understand 
the physical conditions of formation and excitation of H2 and 
the origin of the gas. We suggest two avenues for bringing sig- 
nificant insight into these issues. The first will be to observe 
spectral lines of the CH and CH^ molecules in the optical range 
at high spectral resolution. These molecules are linked to the 
formation and excitation o f H2 (Federman 1982; Mattila 1986; 
ISomerville & Smithlll989l) . The formation of CH is predicted to 
be controlled by gas-phase reactions with H2. The CH is thus a 
good tracer of H2 and their abundances are generally strongly 
correlated. In the presence of a shock, the formation of the CH^ 
molecule through the chemical reaction C^ + H2 needs a tem- 
perature about of 4500 K to occur. Thus, the CH^ molecule is a 
probe of hot and excited media, which could be interpreted for 
our targets as material close to the star, and this would allow us to 
better constrain the excitation of H2. The observation of CH and 
CH^, if present, would provide direct evidence of warm/hot CS 
gas close to the star. Combined with the H2 data, this will allow 
us to explain the chemical mechanisms of formation/destruction 
and excitation of the different molecules in the CS environment 
of HAeBes and to better evaluate the velocity dispersion along 
the line of sight (possibly due to turbulence). 

The second avenue would be to observe pure rotational and 
ro-vibrational transitions of H2 in the infrared wavelength range. 
This would give access to the spatial distribution and mass 
of the observed gas and would help constrain the excitation 
mechanisms. Such infrared observations would also provide in- 
formation about the circumstellar dust, including PAHs. Next- 
generation instrument s Uke VLT/VISIR coul d provide these 
kinds of observations dMartin-Zai'di et al]|2007h . 



All these data would give us the opportunity to obtain a 
global picture of both the structure and the evolution of the CS 
environment of HAeBes. 
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Table 4. Column densities of the different energy levels of H2 when clearly observed in the FUSE spectra of the stars, "g" is the statistical weight of each level. 
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Fig. 8. Overplot of the FUSE spectrum of HD 176386 (thin line) and the synthetic spectrum obtained by combining our photospheric 
model with the results of the Meudon PDR Code (thick line). This plot only presents the portions of the spectrum where the stronger 
absorption lines of H2 ( v = 0, 7 = 0, 1 , 2) are observed. These principal prominent absorption troughs caused by molecular hydrogen 
are indicated throughout the spectrum. Airglow lines are indicated by ffi symbols. 
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Fig. 9. Overplot of the FUSE spectrum of HD 250550 (thin line) and the synthetic spectrum obtained by combining our photospheric 
model with the results of the Meudon PDR Code (thick line). Same legend as for Fig. |8l Here, the photospheric model does not fit 
the data well because of the strong stellar wind lines (for details see iBouret eti 
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Fig. 10. Overplot of the FUSE spectrum of HD 85567 (thin line) and the synthetic spectrum obtained by combining our photospheric 
model with the results of the Meudon PDR Code (thick line). Same legend as for Fig. [8] 
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Fig. 11. Overplot of the FUSE spectrum of HD 259431 (thin Hne) and the synthetic spectrum obtained by combining our photo- 
spheric model with the results of the Meudon PDR Code (thick line). Same legend as for Fig. |8] 
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Fig. 12. Overplot of the FUSE spectrum of HD 38087 (thin line) and the synthetic spectrum obtained by combining our photospheric 
model with the results of the Meudon PDR Code (thick line). Same legend as for Fig. [8] 



